Shot noise of spin polarized electrons 



O 

o 

(N 
X) 

\o 

(N 



A. Lamacraft 

Department of Physics, Princeton University, Princeton, NJ, 08544, USA. Q 
(Dated: September 21, 2003) 

The shot noise of spin polarized electrons is shown to be generically dependent upon spin-flip 
processes. Such a situation represents perhaps the simplest instance where the two-particle character 
of current fluctuations out of equilibrium is explicit, leading to trinomial statistics of charge transfer 
in a single channel model. We calculate the effect of spin-orbit coupling, magnetic impurities, 
and precession in an external magnetic field on the noise in the experimentally relevant cases of 
diffusive wires and lateral semiconductor dots, finding dramatic enhancements of the Fano factor. 
The possibility of using the shot noise to measure the spin-relaxation time in an open mesoscopic 
system is raised. 

PACS numbers: 72.70.-|-m,72.25.-b,72.15.Gd 
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Quantum statistical fluctuations of current and volt- 
age in mesoscopic conductors encode information beyond 
that obtainable from averaged measurements. The sec- 
ond moment of current contains two-particle correlations 
as subtle in principle as the one-particle signatures of 
phase coherence - weak localization, universal conduc- 
tance fluctuations, and so on - revealed by conductance 
measurements. The full counting statistics (FCS)i*^ - 
the probability distribution of charge transfer through a 
conductor - similarly depends on the multiparticle cor- 
relations of a many electron system. From a theoreti- 
cal standpoint the simplest conductor is a scatterer in 
a one dimensional channel of non- interacting fermions. 
Here, the FCS are binomial with a number of attempts 
eVr/h, where V is the voltage bias and t the observation 
time, and a success probability for charge transfer T, the 
transmission coefficient of the scatterer. Thus, although 
the existence of a finite attempt frequency is a signature 
of the correlations imposed by the exclusion principle, 
the 'elementary events' are still described by one-particle 
probabilities. 

It is the purpose of this paper to describe the simplest 
situation in which multiparticle effects enter in a more 
essential way, leading to particular signatures in the sec- 
ond moment and to a novel trinomial statistics of charge 
transfer. We will first describe the physics in an idealized 
conductor, before moving on to the more realistic cases 
of diffusive wires and semiconductor quantum dots. 

We will be concerned with the injection of spin- 
polarized currents into a scatterer where spin-flip pro- 
cesses may occur. The basic phenomenon was identified 
in a recent Letter— for a particular scattering geometry 
and spin-orbit interaction. In fact, the shot noise of spin- 
polarized electrons is generically dependent on spin-fiip 
processes, as these determine the triplet to singlet scat- 
tering amplitude involved in scattering electrons from 
two different channels into the same channel. We pre- 
dict dramatic enhancements of the noise-to-current ra- 
tio (Fano factor) as the injected spin relaxes, raising the 
possibility of using the noise to measure spin-relaxation 
rates. It is remarkable that a charge transport measure- 
ment generically contains information on the dynamics 
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FIG. 1: One-dimensional model illustrating the effect of spin- 
precession on two-particle scattering. Spin up electrons enter 
from left and right, scattering from a potential scatterer while 
precessing in an external magnetic field. 



of spin. The effect of spin flips on shot noise for particu- 
lar spin- valve geometries was also discussed in two recent 
preprints^ 

Consider first the two-particle scattering problem il- 
lustrated in Fig. ^ Spin polarized electrons enter a one- 
dimensional channel of length L from left and right. The 
Hamiltonian is given by _ff = -t- {ez/2)'b ■ a, where 
the unit vector b gives the direction of the external mag- 
netic field, and ez is the Zeeman splitting. In the channel 
at distance £ from one end is a potential scatterer with 
S-matrix 
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with TLs indicating the identity in the spin sector. The 
overall S-matrix of the structure, linearizing the disper- 
sion near the fermi energy, is thus 
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where wj? is the Fermi velocity. The angle of precession 
is determined by the length of time the electrons spend 
in the channel. There are three final outcomes for the 
location of the two electrons: both leave to the left, both 
to the right, or one in each direction. Accounting for 
fermi statistics, and assuming that the electrons have the 
same energy, these probabilities are. 

Both to left or right: TR\Ui^{L - 2£)\^ (1) 
One in each direction: + + 2TR\U^^{L - 2i)f , 

where T = |t|2 = \t'\'^ , R = 1 - T . We see that only 
through the spin rotation are both electrons allowed to 
leave to the same contact. They must form a singlet: 
TR\Uii{L - 2^)p is the probability for this event. If 
the scatterer lies precisely in the middle of the chan- 
nel (L = 2i), reflected and transmitted electrons precess 
through the same angle, so the probability to form a sin- 
glet is zero. Note that assuming that electrons entering 
the channel are polarized in a given direction is equiva- 
lent to having an inhomogcnous magnetic field, vanishing 
outside the channel. A uniform field conserves total spin 
so that there is no triplet to singlet amplitude. Scatter- 
ing due to band distortions in the inhomogcnous field can 
be neglected if we assume the scale of variation is much 
longer than the Fermi wavelength. 

The implication of this consideration for the statis- 
tics of transmitted charge is the following. Evidently no 
charge is transmitted on average, but for |?7|-f (L — 2^)p ^ 
0, there will be uncertainty in the amount of charge that 
has passed through the system, as there is a non-zero 
probability to transport two charges in either direction. 
To make this connection more formal, let us compute the 
FCS for this problem. We obtain the generating function 
x{X,t) = E„^(";T)e*^" of the probabilities P(n,T) to 
pass n charges in time r using the result"'^ 

ln[x(A, r)]^^J dE tr In (l + h{E){SzlS^ - 1)) 



where n{E) — diag(ni, |, nz, j^, n/j |, n^^ |) is the distribu- 
tion function of the incoming electrons. The injection 
of an excess of spin up electrons into both ends of the 
channel corresponds to nL.R].i{E) — 0{fJ-\.i — E) at zero 
temperature, with A/i = — > 0. We find 

x(A,r) = (r2 + i?2 + 2Ti?|[/TT(L-2£)|2 

+2TR\Ui^{L-2(.)\^cos2\)^'""'^ . 

This is the generating function of a trinomial distribu- 
tion with number of attempts A/ir//i and probabilities 
to transfer ±2 and charges given by Eq. The corre- 
sponding noise is& 

S = lim e2(An2)/r = ^A/iTfllC/i t(L - 2^)p , 
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FIG. 2: Fano factor (defined hy S — Fel) for a diffusive 
wire passing a spin-polarized current, as a function of spin- 
relaxation length Ls. Inset: experimental geometry 



which displays a dependence on the spin flip probabil- 
ity |?7||p that mirrors the result of Ref. in a different 
geometry, where the spin rotation was due to spin-orbit 
coupling. Generally, a non-trivial dependence of the noise 
on parameters governing spin-flip scattering may be ex- 
pected when there are two or more channels of incoming 
spin polarized electrons, for then the noise will contain a 
contribution corresponding to two electrons in a triplet 
state and different incoming channels passing to a sin- 
glet in the same outgoing channel. Spin flips may be 
caused by magnetic impurities, spin-orbit scatterers, or 
precession in an external magnetic field if contributing 
trajectories are of different lengths, as in the above ex- 
ample. In the remainder of this paper we will discuss 
the implications of this observation for noise in diffusive 
wires and semiconductor quantum dots, where the cor- 
responding measurements in the unpolarized case have 
been madc^-^. 

In such mesoscopic scatterers, an average over scatter- 
ing matrices is called for. Note that Lesovik's® formula 
S — x^/^Sn^n(l ~ Tn), expressing the noise in terms 
of the eigenvalues r„ of the transmission matrix tH^ does 
not apply in the present situation: the density matrix is 
non-trivial in the incoming channels, so that a rotation 
to the block diagonal S-matrix is not possible. Thus one 
route to the calculation of the averaged noise - averaging 
over the distribution of T,/" - is not available. 

Diffusive wires. We consider the geometry depicted in 
the inset to Fig [3 a diffusive wire of length L and cross- 
sectional area with diffusion constant D. A micro- 
scopic calculation of the shot noise of a diffusive wire was 
presented in Refjllj. Here we pursue a formally equivalent 
approach based on Nazarov's circuit theorjii^. Current 
correlators may be obtained from the Keldysh Green's 
function satisfying 



- - ra — J, ) GA(r, r't, t') = 5{r - r')S{t - t') 
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where denote the Pauli matrices in Keldysh space, and 
Jx = —ie{dx — 'dx)l'2.m is the x-component of the one- 
particle current density operator . With this definition 
we can write the generating functior-'"^ 



x(A,r) = (Tif exp ( ^ / dt j dv^T^Jx^! 



exp I -- y dA'tr{T3J^GA'} j , 



where the trace is understood to be over the Keldysh and 
spin spaces, as well as space and time indices. Tr- denotes 
time ordering along the Keldysh contour. In the above we 
depart slightly from the usual formulation in that we take 
the 'counting field' A to be constant throughout the con- 
ductor instead of defining a surface at which we measure 
the current. This does not affect the result as the zero- 
frequency correlator {I(x)I(x')) is independent of a; and 
x' by current conservation. Using standard quasiclassical 
methods in the diffusion approximation we arrive at the 
following formulation of the disorder averaged problem 
(after Keldysh rotation 



dis 



ievDr 



dE ti-'lngxEdxgxE} 



iez /2[h ■ CT, gxE] 
7;—[<Jig\Ecri,9\E] 

ZTs 



0, 



(2) 



where D is the diffusion constant, v is the density of 
states at the Fermi energy, g{r) — iG{r, v)/ttv is the qua- 
siclassical Green's function, and V = V — iA/2L[r3, ■ • •]. 
tr' denotes a trace over the Keldysh and spin spaces only. 
We have introduced the total spin relaxation rate 1/ts 
due to spin-orbit and magnetic impurities. To find the 
noise we require a solution of Eq. |2] to first order in A 
only. The zeroth order solution has the usual form 



-1 



where Te is a 2 x 2 matrix in spin space, related to the 
distribution function Me by Te = 1 — ^Me- We chose 
the simplest model for spin injection from the left: a 
half-metallic ferromagnet connected to the wire through 
a perfect interface, providing a reservoir of spin up elec- 
trons only at chemical potential /i| . At the right we have 



a normal reservoir with no spin polarization and chemical 
potential ^n- Since we will derive a general expression 
for the noise in terms of Te, extending the calculations 
to more realistic injection scenarios involving incomplete 
polarisation and interface resistances is a matter of ki- 
netics. With the decomposition J^e = Feq^ + - a, 
the corresponding boundary conditions are 

Feq{~L/2) = 1/2 [Feq,„(S - Mt) + i^oqm(S - Mw)] 
F£o(L/2)-Feq,n(S-/iAr) 

Fb(-L/2) - 1/2 [Feq,n(S - Mt) - ^^cqm(^ " Miv)] § 

FEiL/2) = , (3) 
where Fcqm{E) = tanh(i?/2A:T) corresponds to the Fermi 
distribution. The singlet and triplet parts of the distri- 
bution function satisfy 

DV^Feo = 

9 - 2 
DV^Fe - fzh xFe Fe = . 

Ts 

The solution is simplest when the Zeeman term is absent. 
FEoi^) = (1/2 - x/L)FEoi-L/2) + (1/2 + x/L)FEoiL/2) 



FEix) 



Fe{-L/2) 

1 - e"2L/i. 



where Ls — \/ Dts/2. The zeroth order result is thus 

tr' {TiCjEdxgE}^'^^ = ^^"^ involves only the singlet 

part: the average current is unaffected by spin polariza- 
tion. Substituting into Eq. Ogives the first order correc- 
tion 



tr'{Tig£a,5£}(i) = 4a(i^|o + F| - 1) 



(4) 



The constant of integration in Eq. 0] is fixed by the 
fluctuation-dissipation theorem in equilibrium. We thus 
obtain ln[x(A)]/r = iXI-\^S/2+- ■ •, where / = GA/i/2e 
is the average current (G = 2e^vDA/L is the conduc- 
tance including spin degeneracy) and the noise is 



S^^j_^^dxjdE[l- Fl,{x) FUx)] , (5) 



G /-^/^ 

-L/2 

Performing the x and E integrals in Eq. |S1 gives 



S^G 



4kT + A/Lt coth 



/Am 
\2kT 



{Lsim 



2(1 



I 



kT 



(Am/2) coth 



One may verify that this result satisfies the following lim- 



its. As A/i = /i| — /i^v 0, the second term vanishes 
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(originating from the triplet part) and we are left with 
S — 2GkT, as required by the fluctuation-dissipation 
theorem. As ^ oo and at T = 0, the noise is 
S — (l/3)e/. This is the familiar result for unpolarized 
electrons^", indicating that without spin relaxation there 
is no signature of spin polarization in the noise. As Ls 
decreases from infinity to zero, the Fano factor - defined 
hy S = Fel - changes from 1/3 to 2/3 (see Fig.|21). Note 
that although we used a microscopic quantum mechani- 
cal approach, the apperarance in Eq.jSjof only quantities 
from kinetic theory indicates that phase coherence is not 
necessary to observe these effects^^ 

Quantum dots. Spin-dependent transport phenomena 
have recently been the subject of intense investigation 
in lateral semiconductor quantum dotsi^^iStii^. Again, 
relaxation of injected spin-polarized current may be due 
to spin-orbit coupling or precession in an external field. 
As in the case of the diffusive wire, we will consider the 
shot noise in a two terminal set-up where a spin-polarized 
current is driven through one of the terminals. In the 
spin degenerate case, we have the well-known result F = 
NlNr/N'^, valid when N = Nl + Nr :^ 1, where Nl 
and Nr are the number of fully open channels in the 
left and right leads^iiSiSS. Working in the limit of large 
channel number allows us to ignore the effects of weak 
localization and mesoscopic fluctuations. 

We will find the noise in this case by a direct average 
of the expressionSiiSS 



Ah 



J dEtY{{S''AS-A){l+J'E){S^AS-A){l-TE)} 

(6) 

A = (NlPl - NR'Pn)/N^, where Vl/b. are projectors 
onto the channels in the two terminals. The trace is 
over the orbital channels in the leads, as well as spin. 
An average of Eq. |H1 in the N 1 limit using standard 
techniques yields^S. 



S = ^ J2 [ dE{V-{l + Teo.))m(V ■ (1 - Tep))m 

+1^ / dEtY{A{l + J^E)Ht-J^E)}- 

This is the analog of the result ||SJ) for the diffusive wire. 
Gd = {2e'^ /h)NLNii/N is the conductance of the dot. 
Greek letters index channels in the leads. Vlm is the 
zero-dimensional diffuson given by Poo = {^''^[^apSjj^]), 



T>LM = (tr[(TL5a/3crA/5^„]) L,M = 1,2,3, where the 
trace is over the spin indices. The matrices H ± J- are 
understood to be resolved into the singlet-triplet basis as 
before prior to multiplication by the diffuson. We will 
consider two cases: the effect of an external magnetic 
field only (no spin-orbit coupling) , and the effect of spin- 
orbit coupling only. Following the notation of Ref. 0, 
the two cases correspond to. 



N + iezS-h 



In the above {Sk)lm = ~i^KLM are spin-1 operators. 
cf' is the rate of spin relaxation in the x — y plane, 
whereas e^^ governs the relaxation of the z-component 
(both these and ez now expressed in units of A/27r?i, 
where A is the single-particle level-spacing of the closed 
dot). Straightforward calculation using the distribution 
function ^ gives the T = Fano factor 



NlNr , 1 



iV2 



-Nf 



1 



1 
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in the first case, where 6 is the angle between s and b. 
For the second case, we have 



NlNr 



1 



-Nf 



iV2 '^^ (AT + 2eSO)2 



1 



— sm 



{N 



^so 



+ eSO)2 



where 9^ is the angle to the z-axis. A feature of this 
two-dimensional system is that, since e^*^ ^ the de- 
viation of the Fano factor from NlNr/N'^ depends on 
the direction of polarization. In all cases the Fano factor 
of a symmetrical dot {Nl — Nr) increases from 1/4 to a 
maximum value of 3/8 with increasing spin-flip rate. 



Acknowledgments 

This research is supported in part by the David and 
Lucille Packard foundation, and by the National Science 
Foundation under Grant No. PHY99-07949. 



* Electronic address: |alamacra@princeton.edu' 

^ L. S. Levitov and G. B. Lesovik, JETP Lett., 58, 230 

(1993); L. S. Levitov, H. Lee, and G. B. Lesovik, J. Math. 

Phys., 37, 4845, (1996); L. S. Levitov, cond-mat/0210284 
^ Quantum Noise in Mesoscopic Physics, edited by Y. V. 

Nazarov (Kluwer, Dordrecht, 2003). 
^ J. C. Egues, G. Burkard, and D. Loss, Phys. Rev. Lett. 



89, 176401 (2002).^ 

" E. G. Mishchenko, | cond-mat/0305003| 

^ W. Belzig and M. Zareyan, cond-mat /0307070" 

® The definition of the noise used here - a factor of 2 smaller 

than the convention used in e.g. Ref. 21 - is the natural 

one in the context of PCS. 

M. Henny et al., Phys. Rev. B 59, 2871 (1999) 



5 



* S. Obcrholzer et al., Phys. Rev. Lett. 86, 2114 (2001). 
^ G. B. Lcsovik, JETP Lett. 49, 592 (1989); M. Biittiker, 

Phys. Rev. Lett. 65, 2901, (1990). 
^° C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 
1889 (1992). 

" B. L. Altshuler, L. S. Levitov, and A. Yu. Yakovets, JETP 
Lett. 59, 857 (1994). 

Yu. V. Nazarov, Ann. Phys. (Leipzig) 8, 193, (1999). 
A. L Larkin and Yu. N. Ovchinnikov in Noneqmlibrium 
Superconductivity, edited by D. N. Langenberg and A. L 
Larkin, North-Holland, 1986. Appendix to Chapter 11. 

" K. E. Nagaev, Phys. Lett. A, 169, 103 (1992). 

1^ L L. Aleiner and V. L Fal'ko, Phys. Rev. Lett. 87, 256801 
(2001) 



^® P. W. Brouwcr, ,L N. H. J. Cremers, and B. L Halperin, 
Phys. Rev. B 65, 081302 (2002) 

D. M. Zumbuhl, J. B. Miller, C. M. Marcus, K. Campman, 
and A. C. Gossard, Phys. Rev. Lett. 89, 276803 (2002). 
^® J. B. Miller et al, Phys. Rev. Lett. 90, 076807 (2003) 
R. A. Jalabert, J. L. Pichard, and C. W. J. Beenakker, 
Europhys. Lett. 27, 255 (1994). 

Yu. V. Nazarov, in Quantum dynamics of suhmicron struc- 
tures, edited by H. A. Cerdeira, B. Kramer, and G. Schon, 
NATO ASl Series E291 (Kluwer, Dordrecht, 1995). 
2^ Ya. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 
M. L. Pohanski, M. G. Vavilov, and P. W. Brouwer, Phys. 
Rev. B 65, 245314 (2002). 



